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Abstract. We demonstrate the generation of classical incoherent light with electronic 
control over its temporal characteristics and photon number distribution. The 
tunability of the temporal coherence is shown, under both classical and quantum 
detection, through second order correlation (G'^(t)) measurements. The tailoring 
of desired classical photon number distributions is illustrated by creating two 
representative light sources - one thermal and the other a specific classical, non- 
Gaussian state. Such generation of classical light sources, quite different from 
existing natural light sources, is likely to be a useful resource in quantum information 
processing. As a particular application in this direction we outline how a non- 
Gaussian state generated in this manner may be mixed with an appropriate non- 
classical Gaussian state at a beamsplitter, to generate non-Gaussian entanglement. 
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1. Introduction 

Understanding light-matter interactions and mechanisms of hght generation are of 
importance in elucidating many physical processes. For most processes, a combination 
of a statistical approach along with Maxwell's classical theory of electromagnetism 
suffices. However, several optical processes demand a quantum mechanical approach 
to the electromagnetic field [1, 2, 3]. The quantum theory of hght makes a clear 
distinction between different sources of radiation such as coherent, thermal, and single 
photon sources [4] . While the radiation from the former two sources may be explained 
classically, the light from the latter can only be explained quantum mechanically. In 
other words, the former two states are deemed classical, and the latter non-classical. 
Due to their quantum nature, non-classical sources of light display several counter- 
intuitive features evoking considerable interest in them. Such sources of photons have 
been generated experimentally by several techniques [5] and their possible applications 
in different fields, especially quantum information, are well explored [6]. However, recent 
years have witnessed an emergence of interest in classical sources of light such as thermal 
light. These have been used in intriguing applications like ghost imaging [7]- [10], and 
interfcromctry based experiments as in Refs. [11, 12], and also used as an ingredient in 
the creation of non-classical states of light [13, 14, 15]. 

In this paper we demonstrate and elaborate upon a method of creating classical 
incoherent light quite distinct from that obtained from existing natural light sources. 
The method utihzes the fact that an acousto-optic modulator (AOM) may be effectively 
used to introduce phase and intensity fluctuations to light rapidly and accurately [16, 17]. 
Starting from coherent laser light, we generate incoherent light with desired coherence 
times and desired intensity statistics, or photon-number distributions. The motivations 
for this are several. Not only does it offer an alternative to the current method of using 
a rotating ground glass plate [18, 19] to generate pseudo-thermal light, it offers a robust 
and ffexible procedure for producing tailored classical light sources with predetermined 
photon emission rates, which can be used in experiments such as photon addition and 
subtraction [13]-[15] ,[20]-[26]. Several interesting applications now become possible. For 
example, as we show in this paper, classical non-Gaussian states with tailored photon 
number distributions (PND) created by this technique may, in turn, be used to produce 
states with tailored non-Gaussian entanglement. This is of importance in the quantum 
information theoretic context, where recent findings suggest non-Gaussian entanglement 
to be advantageous over Gaussian entanglement [27, 28, 29]. 

This paper is arranged as follows. The experimental setup used in this study is 
described in Section 2. Section 3 describes the methodology for the generation of classical 
incoherent light with tunable temporal coherence. This is achieved by imparting random 
phase shifts to coherent light through the acousto-optic interaction, in a Mach-Zehnder 
interferometer (MZI). Measurements of the second-order correlation function, G^(t), 
both under classical (intensity-intensity correlation) and quantum (photon coincidence) 
detection, show the tunability of temporal characteristics, and the equivalence of the 
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two forms of detection. In Section 4, we generate classical light sources with the 
desired PNDs introducing intensity fluctuations to light by suitably modulating the 
diffraction efficiency of the acousto-optic modulator (AOM). Two incoherent states, 
one thermal and the other a classical non-Gaussian state were created as illustrative 
examples. In Section 5 we discuss the possible use of such tailored classical non-Gaussian 
state in producing tailored non-Gaussian entanglement. Section 6 summarizes the work 
presented in this paper. 

2. Experimental Setup 




The schematic of the experiment is given in Fig. 1. A Mach-Zehnder interferometer 
(MZI) was used, with an AOM in each of the two arms, arranged such that the first 
order diffracted light proceeded further on, while the undiffracted light was blocked. A 
coherent light beam at 767 nm from a diode laser (Toptica), the intensity of which could 
be adjusted using the variable neutral density filter (VNDF), was amplitude divided at 
the 50:50 beam splitter BSl. After traversing the two arms of the interferometer, the 
beams were combined at BS2. The two AO Ms were driven by individual Versatile 
Function Generators (VFG, Toptica) which operated at 80 MHz radio frequency (rf); 
they were both referenced to a common stable clock. Fine electronic control over the 
phase and intensity of the light was achieved by programming the phase and amplitude 
of the radio-frequency output of the VFG using a Lab VIEW interface. One could tailor 
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any distribution of phase and intensity fluctuations (on time scales ranging from few 
hundreds of nanoseconds to seconds) in the rf electrical signal being fed to the AOM. 
These fluctuations are transferred to the diffracted light by the acousto-optic interaction. 
The diffracted light beams of the two arms of the MZI interfere at BS2. Light emerging 
from one of the exit ports of BS2 was directed into an input of a 50:50 flber-splitter, 
and then onto two detectors, Dl and D2. Measurements at the two detectors were used 
to determine the second order correlation (intensity-intensity correlation) function. For 
the case of classical detection of light, the laser was operated at a power of around 
30 mW and Dl and D2 were two fast photodiodes (Thorlabs PDIOA-EC). For the case 
of quantum detection, the laser light was strongly attenuated, photodiodes Dl and 
D2 were avalanche photo-diode (APD) based single photon counting modules ( SPCM- 
AQR-15 Perkin Elmer) where an incident photon results in a TTL pulse with a detection 
efficiency of 65 % at 767 nm. The outputs of Dl and D2 were stored in a PC using data 
acquisition systems. In the case of classical detection, a digital storage oscilloscope was 
employed, while for quantum detection, two counters on a data acquisition card (NI 
M-series PCI-6259) were used. 



3. Generation of Incoherent Light by only phase modulation in an MZI 
setup 

In this Section, we describe the creation of incoherent sources of light by imparting phase 
jumps to light by means of the AOMs in the MZI (Fig.l). The use of a Mach-Zehnder 
interferometer was motivated by the facts that, i) phase changes imparted to light can 
be discerned only in an interferometric setup, and ii) it enables the creation of a source 
with intensity fluctuations even though only phase jumps are imparted to light. 

3.0.1. Theory The action of the MZI in Fig.l may be mathematically represented by 
the transformation matrix M = B.^.B where B ( the operation of a beam splitter) and 
$ (the action of the two AOMs) are given by 

B = ^ ( ' ' \ and 4 = ( ° „ \ (1) 

^2 y -1 1 y y e^'^^W j ^ ' 

Here 0i(t) and 02(^) are the phase shifts imparted to light at AOMs 1 and 2, respectively. 
Thus, in terms of the scalar wave field, i^i, entering one of the ports of BS\ (and with 
no input at its other port) the output at BS2 is given by 

(2) 

where 1,2 represent the input ports of BSl and 1', 2' the output ports of BS2. It is clear 
that the second order intensity self-correlation at the output ports is given by 

^ {EmE,{t)),{Et{t + r)E,{t + T)), ^ ^ 
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where i can take values l',2'. In terms of field operators Gf^{T) is given by 



A^t + r)Mt + r) (V'l Ai^t)A{t) 



where 





(4) 



(5) 



with 1^) being the input state, Oi, 02 the annihilation operators of the input modes, 
and Aii,A2' those of the output modes, and M is as defined earlier. In our setup, 
IV') = \Ei) (8) |0), where is a coherent state. Simplification of Eqs. 3 and 4, for the 
case of complete phase noise, (i.e., uniformly distributed noise, with phase spanning the 
entire circle) both yield 

GUr) = 1 + {cos{5m cos(50(i + r))), (6) 

where 6(f) — 4>i{t) — 4>2{t)- Note that we have adjusted the MZI to match the dynamical 
paths in both arms and, therefore, the phase difference is the difference in the phases 
imparted at the two AOMs. 

The temporal coherence of light, emerging out of the port 2', is determined by 
intensity correlation technique in a standard correlation setup developed by Hanbury- 
Brown and Twiss [30]. Eq. (6) can be easily cast into the form 

GI,2.{t)^1 + 0.5^{t), (7) 

where ^(r) represents the probability of temporal overlap of S(l){t) and + r). In 
general, for complete phase noise given to AOM 1 and 2 and phase jumps occurring at 
time intervals t given by independent distributions Pi{t) and P2{t) respectively, we have 

^ J-it-r)md{t) J-{t-r)P2mt) 

jrtp,mt) ■ rtp^mt) 

For the special case where one of the phases, say 02 (^), is held constant Eq. 8 reduces 
to 

, ^ ir(t-r)m)4t) 

3.0.2. Experimental Results We now present our experimental results. Initially, 
random phase jumps were imparted to AOMl only and at regular intervals, that is, 
Pi{t) — 5{t — Tc). Prom Eq. 9 and Eq. 7, one obtains 

r 

= 1 for T > Tc (10) 

Thus, on imparting phase noise one expects photon-bunching with a zero-delay second- 
order correlation value 1.5, that falls to 1 for long delays. In the experiment, at constant 
time intervals, Tc, the rf electrical signal to AOMl was given a random phase jump. 



G'^,2,(t) = 1 + 0.5 ( 1 -- 1 for T<Tc 
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Figure 2. G'^{t) as function of r for (a) constant dwell time, noise in one AOM, and 
classical detection; (b) exponential distribution of dwell times, noise in both AOMs, 
and classical detection; (c) same as (a), but quantum detection; (d) same as (b), but 
quantum detection. 

distributed uniformly in the interval (— 7r,7r). The experimentally determined values of 
G'^{t) are shown as discrete points in Fig. 2(a) for three different values of coherence 
time, Tc- The continuous curves are G'^{t) as obtained from Eq. 10. 

Next, AOMs 1 and 2 were driven with rf voltages with independent, random phase 
jumps uniformly distributed in the interval (— vr,??), and at time intervals that had 
independent exponentially falling distributions, P{t) = [exp(— t/rc)]/rc, with 1ms 
< t < 100 ms and with a mean, ~ r^. The results for this case are presented in Fig. 
2(b), for three values of r^. As in the previous case, the agreement between theory and 
experiment is very good. 

The above measurements, which utilized classical detection of intensities with 
photodetectors, were repeated at single photon detection levels. The laser light was 
strongly attenuated, so that, on an average, every 30 /xs, there was a 10% probability 
of detecting a photon. Thus, there was, on an average, less than one photon in the 
interferometer at any instant of time. The output pulses from Dl and D2 were fed to 
two counters operating at a time bin of 30 /is, that was shorter than the mean interval 
between phase jumps that ranged from 1ms to 100 ms. The values of G'^(r) obtained 
from the experiments for the two cases (phase jumps to one AOM only at constant 
intervals and phase jumps to both AOMs with exponential distribution of interval 
between jumps) are shown in Figs. 2(c) and (d), along with theoretically expected 
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values. The agreement is fair. In order to obtain sufficient statistics at single photon 
detection level, data had to be acquired over long durations. Mechanical instability of 
the interferometer and dark counts are believed to have contributed to the deviation 
from the theoretical curves in the quantum detection case. 

From the above, it is amply clear that classical light sources with the desired 
temporal coherence characteristics may be created. The results also confirm the 
equivalence between classical and quantum detection for classical states of light. In the 
present experiment, we have imparted fiuctuations on the milliseconds time scales due to 
restriction of the data acquisition card. The setup allows for a temporal variation from 
50ns (currently limited to 500ns by USB 2.0 communication) to few seconds depending 
upon the stability of the interferometer. This provides an easy way of creating bunched 
light with long coherence times. Further the bunching can be enhanced and higher values 
of G2'2'(0) can be obtained by using engineered partial phase noise in this interferometric 
setup [17]. 

4. Tailoring Photon Number Distribution by Intensity Modulation of Light 

In this section, we demonstrate the use of a single AOM in the generation of sources 
of classical incoherent light with desired photon number distributions (PND) starting 
from an input coherent state. The diffraction efficiency of the AOM is modulated by 
the addition of amplitude noise with a chosen probability distribution to the input rf 
electrical signal. This, in effect, modulates the transmittivity of a coherent state through 
the AOM, thereby providing a source of classical light with a tailored photon number 
distribution. 

4-0.3. Theory In modulating the diffraction efficiency of the AOM, or equivalently, the 
transmittivity of the coherent state, we attempt to realize radiation state 

r\ao\ °° 



where p{n) is the photon number distribution and \a) is the diffracted coherent state 
whose amplitude varies in time according to T^dctl), \ao) is the diffracted coherent state 
at maximal transmittivity, and {\n)} is the Fock basis. The ensemble in Eq. (11) is 
practically realized by modulating the transmittivity of the input coherent state over a 
sufficient amount of time. The upper limit of p can be taken as infinity if |q;o| is chosen 
to be much larger than the mean \a\. 

We experimentally generate light sources with two specific PNDs as examples. The 
first is the thermal state with an average of n photons, namely, pth{n) [2] with 



P 




n=0 





r{\a\) = VthiH) = ^exp(-|a;|Vn), 



(12) 
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for which the probabihties 

Pth(ri) = (1 - e-^)e-^" with e"^ = n/(n + 1), (13) 
define the PND. The second is the state corresponding to 

V{\a\) ^ V^{\a\) = C>|'expK|a|'], (14) 
with an average number of photons n and a PND specified by the probabihties 

pAn) = — . with C = -• 15 

Note that the state is manifestly non-Gaussian. 

The temporal coherence characteristics of the incoherent light with tailored PNDs, 
as generated above, is described by the second-order correlation function G^(r), given 
by 

G\T) = l + {G\0)~mr) (16) 

where ^(r) is as given in Eq. (9). ^^(o) = 2 for Pth(|a|), ^^(0) = 1.5 for P^d"!) and 
G^(0) = 1 for a coherent state. 




Figure 3. (a) Diffracted liglit PND, in absence of modulation (b) PND vs n for 
different A (c) PND vs n for different C 
Insets are corresponding G^(t). 



4.0.4- Experimental Results For this part of the experiment AOMl of Fig.l was 
switched off and BS2 was removed. Thus, only the light emerging from A0M2 could 
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reach the detectors, which in this case, were APDs. The input laser beam was attenuated 
using neutral density filters such that there were, at most, 90000 counts per second. The 
intensity of the +1** order diffracted light was modulated by feeding an rf electrical signal 
with random variations in amplitude that followed the desired probability distribution 
function. The fluctuations in amplitude were made to occur at random time intervals 
in the range 1 to 100 ms with a mean time of ~ 10 ms, with the distribution of time 
intervals falling exponentially. The counter was set to a time bin of SO/is, and the 
sequence of counts was recorded for a period of 30 minutes to obtain good statistics. 
The number distributions of photons was then determined from the data using a time 
bin of 450/is. 

Initially, no noise was added to the rf signal to the AOM. The emergent light was 
the maximally transmitted (diffracted) coherent state \ao); its PND, obtained from 
our measurements was as shown in Fig. 3(a). It gave a Poissonian distribution with a 
mean photon number of 42. The second-order correlation function, G^(t) for this light is 
shown in the inset. When RF power fed to the AOM was varied in the form 7^th(|Q;|), the 
emergent hght was found to have the PND as in Fig. 3(b), and a second-order correlation 
function as shown in the inset. The photons show bunching, and the zero-delay second- 
order correlation has a value 2, seen for thermal light. On similar lines, when RF power 
fed to the AOM was varied in the form P^(|q;|) , the PND obtained for the emergent 
light was as shown in Fig. 3(c), and the correlation function as shown in the inset. In 
all cases, the theoretically expected quantities are shown by the continuous curves. The 
good agreement between theory and experiment underlines the reliability and efficacy of 
this method in generating tailored classical light sources with desired PND and temporal 
coherence characteristics. Generation of classical non-Gaussian states with desired PND 
in a predetermined manner promises to be useful in various contexts, as discussed below. 

5. Discussion 

The introduction of phase and intensity ffuctuations to coherent light using a single 
AOM gives rise to incoherent light with a clean spatial distribution, and permits the 
electronic tuning of temporal characteristics. As light with clean spatial profile with 
phase and intensity fluctuations is a requirement for experiments on single-photon 
addition/subtraction to thermal light [14], this AOM-based technique is likely to 
flnd immediate application here. Further, the tailored classical non-Gaussian states 
generated by the techniques discussed so far could be used in conjunction with photon 
addition/subtraction to create tailored non-classical non-Gaussian states. As the AOM 
offers complete ffexibility of providing phase and/or intensity ffuctuations with different 
desired probability distributions and on different time scales, the device can be used to 
generate classical incoherent light having properties quite different from existing light 
sources, opening up new flelds for exploration. Another feature, likely to prove useful, 
is the ability to control the mean photon number with ease, even during the course of 
a measurement. This has been shown for the two speciflc cases in section 4, where we 
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have changed the parameters A and ^ to alter the mean photon number. Now some 
specific apphcations are outhned below. 

5.1. Generation of non-Gaussian Entanglement 

Here we show how the tailored PNDs of the kind generated above may be used for 
producing novel forms of entanglement. A simple experimental scenario is discussed. 

As is well known, a beamsplitter preserves non-Gaussianity [31], and non- 
classicality [2, 32]. Further it is known to generate entanglement from input non-classical 
states [33]. This feature of the beamsplitter may therefore be utilized to create non- 
Gaussian entangled states [24] of the form 



by appropriately choosing a Gaussian state po and a tailored non-Gaussian state pnG- 
Here U^s is the beamsplitter unitary. 

Consider the situation where a single-mode squeezed state psq enters one port of a 
50:50 beamsplitter and a tailored non-Gaussian state enters the other. Note that 
here we may also tailor the non-classicality of pnc through processes such as photon- 
addition. Let V — diag \ {e'^, e~^) denote the variance matrix of the squeezed state with 
squeezing parameter /i, and the tailored non-Gaussian state with average photon 
number n. The resultant state at the output of the beamsphtter is definitely entangled 
when 2n + 1 < [34] , while also remaining non-Gaussian. Clearly, this a realization 
of a non-Gaussian entangled state. The method is effective even if the the initial non- 
classicality were to purely reside in p^Q. For instance, say p = {pq is the ground state), 
and pnG is a photon-added tailored non-Gaussian PND, the resulting p"^ is definitely 
non-Gaussian entangled [35]. In other words, we may choose the initial non-classicality 
to reside either in pc or in p^Q, or in both of them, while simultaneously tuning the 
non-Gaussianity of pnG, to generate the desired non-Gaussian entanglement. 

6. Conclusion 

To conclude, we have experimentally demonstrated a method for the creation of 
incoherent classical light using AOM where the temporal characteristics, coherence time, 
photon number distribution function and the mean photon number is electronically 
tuned. Specifically, creation of pseudo-thermal light and non-Gaussian classical light 
have been chosen as illustrative examples. Possible apphcations of such tailored light 
sources have been discussed, like generation of tailored non-Gaussian entanglement. The 
present proof-of-principle experiments, which display fluctuations on the milliseconds 
timescales, may be easily augmented with currently available technology to higher speeds 
(nanosecond timescales) . 




(17) 
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